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SOLVING  THE  BOUNDARY  LAYER  PROBLEM 
V.  Ya.  Shkadov 

The  movement  of  an  incompressible  viscous  fl&ld  in  a  boundary 
layer  is  described  in  Kotschin  °t  al.  [1]  by  the  equation  for  the 
stream  function  ^  (x,y) 

d.|)  a*  1*2.  a  Hail 

0; /  bx  by  bx  by*  dx  V  by* 

with  the  boundary  conditions 

d\b  t)ih  /  \ 

*  =  0.  ^-  =  0  !/  =  0,  —  -!/*•»■*«  '  (2) 

If  the  velocity  of  the  fluid  U(x)  on  the  external  boundary  of 
the  boundary  layer  is  given  in  the  form  of  a  power  series  U(x)  =  aD  + 

+  ajx  +  a2x*  +  . ..,  the  solution  for  if/  (x,y)  may  also  be  found  in  the 
form  of  a  series  with  respect  to  x,  the  coefficients  of  which  must  be 
found  by  numerical  integration. 

At  present  Curie  [2]  has  carried  calculations  out  to  the  x11  term 
for  bodies  with  blunted  fronts  and  symmetrical  with  respect  to  the 
onflowing  stream.  When  using  this  solution  one  must  limit  himself  to 
the  first  few  terms  of  the  U(x)  expansion,  but  this  is  not  always 
sufficient,  especially  for  asymmetrical  bodies.  In  these  cases  one 
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must  use  a  solution  which  cannot  represent  U(x)  as  a  rapidly  converging 
series. 

The  author,  himself  [J]  has  pointed  out  the  possibility  of  pro¬ 
ducing  a  solution  to  the  boundary  layer  equations  definable  by  the 
dimensionless  combinations  U«x/U,  U»  «x2/U,  U"  'xs/U,  ... ,  In  which 
the  primes  Indicate  differentiation  with  respect  to  x. 

We  will  Introduce  a  new  Independent  variable  q  =  q  (x,y)  and  the 
unknown  function  f  (x,q)  so  that 


(x<  y)  —  V  u  vx 


'  (*•*»)  (’i  =  *V^) 


For  f  (x,q)  Is  derived  the  equation 

<?i  t  !)■!  .  vx\  t<yy  ,  j_  *>sn_  ( df  d'f  _  iL  t IL\ 
- 5i? + T t oif  *r  u  V~\dnl  •  2  •  bx  d1?' 


O) 


with  the  bounding  conditions 


7-o, 


Of 

0'\ 


=  0 


n-o. 


df  . 


(4) 


Let  us  examine  the  equation 


gy  ‘  i  td-/  u'x  r 
00 

=  ^  i'Pi  —  PlPi^PiA 
i== I 


V, 


:-i'  \  e>r]  dp{bx\ 


jL<m 

8  Pi  <pt?) 


In  which  the  function  depends  on  q,  pr,  p2,  ... 

We  will  look  for  a  solution  to  this  equation  which  satisfies  the 
boundary  conditions  (4)  .  If,  In  the  space  of  the  variables  pt,  p2, 
...,  function  U(x)  parametrically  gives  a  curve  with  the  relationships 
Pi  =  Urx/U,  p2  =  U* 'x2/U,  It  Is  easy  to  see  that  along  this  curve 

xip-  =*>4  —  PiPi  +  Pi+i  (>  =  1.2 )• 

-  J  •  ax 

therefore  the  right  half  of  Eq.  (5)  represents 


/a/  _£_=/ _ 'oj_  ay_ \ 

\drj  dxdr\  dx  drf) 
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and  Eq.  (5)  coincides  with  (3)  . 

Consequently,  in  this  special  case  function  f(q,  px,  p2,  ...) 
found  by  the  solution  of  Eq.  (5)  with  boundary  conditions  (4)  will 
describe  the  flow  in  the  boundary  layer. 

Two  solutions  for  Eq.  (5)  may  be  produced:  expansion  into  a 
series  with  respect  to  px,  p2,  ...  and  into  one  with  respect  to  px-l, 
pz,  ...  The  first  solution  corresponds  to  boundary  layer  flow  begin- 
ning  at  the  sharp  edge  on  which  px  =  0,  p2  =  0,  ...,  while  the  second 
describe  a  boundary  layer  beginning  at  a  critical  point  at  which 
Pi  ■  1*  P2  =  0 ,  . . . 

The  coefficients  of  these  series  are  functions  of  q  and  are  found 
by  numerical  integration  of  ordinary  differential  equations.  The 
necessary  calculations  ware  done  on  a  "Strela”  computer.  They  show 
that  the  coefficients  of  the  series  rapidly  diminish  as  i^  grows, 
therefore  in  what  follows  the  terms  containing  p4,  p5,  ...  are  con¬ 
sidered  too  small  and  are  dropped.  Por  bodies  with  a  sharp  forward 
edge 

/  =  /oo  T  S/lu/’J  4”  8"  (Aiu/’l2  ~r  fiiPs)  T  S3  (fatiPl*  "T*  hlPlpl~T  f*iP*Y  1+  f5) 


Substituting  £  into  Eq.  (5)  and  setting  expressions  with  various 
combinations  of  px,  p2,  and  p3  equal  to  zero  we  may  derive  equations 
for  f±k. 

Because,  however,  numerical  integration  necessitates  awkward 
right  sides  and  consequent  storage  of  much  information,  flk  was  com¬ 
puted  by  a  different  method.  Let 

U  (x)  —  1  ax  bx 5  ~  c j-3 

It  can  be  found  that 


(a,  b,  e  =  const) 


1  I  1 
:  =  T  “  T p- ' 


0  p'Ji  • 


..  _  JL  ( JL  *  \ 

—  d  \  2  Pi~  2  A/  • 


bx 


,  1  l 


(7) 
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By  expanding  into  a  series  with  respect  to  x  we  will  find  the 
solution  to  Eq.  (j5)  with  the  boundary  conditions  (4)  In  the  form 

/  =  /» -f  («’* i  -r  w=)  Xs  -r  («^-i  —  «Mr.  -f  c*»)  (8) 

The  functions  fQ  (q)  ,  gi(q),  hi(q),  ...  satisfy  the  equations 


/o  Tv  Z0/0"  ~  0 

,  1  .  -  r  3  5 

**  .  2  ^  TT/#'?t  =  -l  t/V5  — y/o// 

*  ”  -  T  *  ~  A>V  -  4  />»  = 

—  -2  ,  t  1 

-  M,‘  “  -  -/•  Hi  —  Y  (/»?«  /#'ei)  t  1  —  /«'=  -r  —  /»/«' 


with  the  boundary  conditions 


A  =  o, 

/»'■-*  i. 

If  we  assume  that  px,  p2,  and  p3  are  small  enough.  It  Is  possible 
to  expand  ax,  bx2,  and  cx2  Into  series  of  the  form  (6)  by  using  (7). 
Substituting  expressions  for  them  into  (8)  ,  collecting  the  terms  with 
the  same  px,  p2,  p3  combinations,  and  comparing  with  (6)  we  finally 
find  For  shearing  stress  at  a  wall  r  =  |i  d2^/dy2  when  y  =  0 

we  have  _i_ 

i\(J *  \  * 

p-1  J  r  =  /'<x>  (0)  -r  8pi/”to  (0)  +  S2  (0)  «f  PiTtk  (0))  + 

-r  S3  (PlYyo  (0)  +  Pi&hi”  (0)  4-  (0))  T  S-*  (0)  -f  (0)  + 

t Pip3/H*2  (0)  -f  (0))  -r  8*  (/>i6/'w  (Q)  t  P*pJ\i  (0)  -r  p£p*fmit  (0)4*  ( 

t  Pipi*i"h3  (0)  *r  PiP*f’  m  (0)) 

where  the  second  derivatives  of  f  *  *  have  the  following  values: 


Ao  (0)  = 

0.33206, 

/■  10  (0)  = 

0.19290, 

im»  (0)  * 

—  6.03129 

A> (0)  = 

—  0.00318, 

/'» (0)  = 

0.012-54, 

/'m  (0)  — 

0.00132 

;*32  (0)  = 

0.00008, 

r«(0)= 

—  0.00679, 

/'« (0)  = 

—  0.00074 

;V.(0)  = 

—  0.00003, 

(0)  = 

—  0.00002, 

/'ic  (0)  = 

0.00423 

/•»  (0)- 

0.00053, 

/'»2 (0) = 

0.00002, 

/'m  (0)  — 

0.00002 

/*S4  (0)  = 

0.000001 

/o'  —  0,  £i  =  0,  »i  —  0,  .  .  .  wh»n  ^  =  o 
A'i'“*0,  ///  — 0,  ...  J),-oo 


(10) 


-4- 


ptd-TT-63-285/1+2+4 


The  calculations  from  formula  (ll)  coincide  well  with  the  results 
shown  In  Gortler  and  Witting  [4]  and  Terrill  [5j.  To  exemplify,  break¬ 
away  of  the  boundary  layer  U  =  1 -x2  takes  place  when  x  =  0.64,  closely 
matching  the  value  of  x  =  0.67  [5] . 

In  order  to  examine  the  boundary  layer  on  bodies  with  blunted 
foreparts  we  will  introduce 


i:'t.  ■ 

IT  ~ 


Vr- 
U  ' 


/U'x  \ 

**  =  —  ~  1  •  x'-  =  ~[T  ~Z\U  ~ 1 J  1 


uv  r 

**  =  "77 - 6 


(12) 


Let  U  =  x  (l  -r  ax2  +  bx4  +  cxe)  ;  then 


i  1 1  i  ,  i  \  .  .  i  f i  *  ,\ 

■=  77  \T*‘  — TXt +  Ic  /  ’  = 

i  / 1  \  i  i  i 

~U\i»  T,j  •  />=  1— v-x.-f-g-x,- 4g*a 


(15) 


Using  the  solutions  found  by  expanding  with  respect  to  x  (Curie 
[2]),  we  will  by  the  above-described  method  obtain 

•  ,  (14) 

f  =  FM  4-  x,Fl0  -r  ~  +  --  • 

and  for  shearing  stress  at  a  wall 


F>-( 


vtf» 


j  t  =  F’m  (0)  xtF'i u  (0)  -f  ritF’o>  (0)  -f  (0)  + 

•f  »|V*  (0)  -+  x^F’z,  (ti)  ~  t3F’,:  (0)  -f  xt*F’ 4o  (0)  +  xfx.F’a  (0)  ■ 
+  rtXjF'r.  (0)  -f  xJF’u  (0)  -f  (0)  -f  (0)  + 

-f  xjVa F’a  (0)  -f  TXxJfa  (0)  J-  (0) 


(15) 


Here 


/-c„(0)=  1.232588, 

(0)  =  — 0.01731, 
0,00111, 
f (0)  =  —  0.00058, 
(0)  =  0.00107. 

F'u  (0)  =  0  00008 


/••„(0)—  0.52W5, 

/•'jo  (0)  =  0  01008, 

r-.o  (0)  =  —  0.00027, 
F’o(0)  =  —  0.00011 , 
F’-.i  (0)  =  0.00010, 


/•'a  (0)  =  —  0.06739 
(0)  =  0.00758 

F’a  (0)  =  —  0.00229 
F"  „  (0)  =  O.O027C 

/•*M(0)=  0.00018 
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By  using  (15)  the  shear  stress  and  point  of  breakaway  are  cal-  , 

i 

culated  for  various  cases  described  In  Curie  [2]  and  Terrill  [6]; 

In  doing  so  good  coincidence  Is  ascertained..  By  way  of  example  we 
cite  below  values  of  the  magnitudes  j 

calculated  for  U  =  U0  (x-x3  + 
also  borrowed  from  Curie  [2], 


■( 


lit 

{/'*  '  oy- 


O.O789  x5)  according  to  formula  (15)  and 


x  .  0.2  0 A  0.5C  o SA 

T  =  1 .  1S9  1  .u27  0.7 15  0.378 

T  1.18 J  ■  1.027  0.712  0.309 


0 . COS  '  0 . coo 

8  _  —  by  fonult  (IV) 

^  frc  Curia  [fj 


Received  January  24,  1962 
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